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In this note, we solve an extended version of the M = 1 super Yang-Mills theory 
with gauge group U(V), an adjoint chiral multiplet and Nf flavors of quarks, by using 
the Af = 1 microscopic formalism based on Nekrasov's sums over colored partitions. 
Our main new result is the computation of the general mesonic operators. We prove 
that the generalized Konishi anomaly equations with flavors are satisfied at the non- 
perturbative level. This yields in particular a microscopic, first principle derivation of 
the matrix model disk diagram contributions that must be included in the Dijkgraaf- 
Vafa approach. 
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1 Introduction 



The recently developed microscopic formalism to M = 1 supersymmetric gauge theory 
[1, 2, 3] is a first-principle approach that allows in principle to solve rigorously a 
general M = 1 gauge theory in the chiral sector. It is based on the direct calculation 
of the relevant path integrals and relies heavily on Nekrasov's instanton technology [4], 
suitably adapted to the Af = 1 context. The goal of the present paper is to develop 
the formalism in the case where fundamental quark flavors are present. Including 
quarks brings several new interesting features that we shall explain in details in the 
following. 

The model we consider is the U(iV) supersymmetric gauge theory with an adjoint 
chiral multiplet X and chiral multiplets Qf and Q* , 1 < / < Nf, in the fundamental 
and anti-fundamental representations respectively. We shall restrict ourselves to the 
cases N t < 2N. 

The lagrangian is given by 

C = Re J d 2 9 ^- Tr W a W a + J d 2 9d 2 9 TrX*e 2V X 

+ j d 2 6d 2 6Q^e 2V Q f + j d 2 6d 2 6& f ~e 2V Qf + 2NRe j d 2 6 W , (1.1) 

where 

is the gauge coupling constant and 

W = 16^ Tr ( V ( X ) WaW «) + + T Q f m f f ~(X)Q f . (1.3) 

The polynomials V(z), W{z) and rrJ ~{z) are parametrized as follows: 

NV(z) = Y,kt k z k -\ (1.4) 

k>2 
fc>0 

m f f (z)=Y, m lj zk - (L6) 

k>0 

It is useful to introduce 

i 

U(z) =det (m(z)) =U Y[(z-b Q ) (1.7) 
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and 

fe>l 

with 

*i = ln(i/og) ■ (1.9) 

The coupling g is the instanton factor in the model. In the asymptotically free case 
(iVf < 2AT) it is expressed in terms of the dynamically generated scale A as 

q = A 2N ~ Nf , (1.10) 

whereas in the case of vanishing (3 function (jV f = 2iV) we have 

q = e 2mr = e^-^lf . (1.11) 

The set of couplings {tfc}fc>2, {fi'fe}, W{ ?} and {b Q } will be denoted by t, g, m and 
6 respectively. In the most standard case, t k = for k > 2, but as explained in [3] 
it is actually very natural and convenient to consider the more general theory with 
arbitrary couplings t. 

The chiral ring of the model is generated by the operators 

TiX k , TrW a W a X k , T Q f 'x k Q f (1.12) 

and our goal is to compute, from first principles, the corresponding expectation values. 
These expectation values are multi-valued analytic functions of the parameters t, g, 
m and q. The multi-valuedness comes from the existence of several distinct vacua in 
the theory. At the classical level, the most general supersymmetric vacuum, which 
is obtained by extremizing the tree-level superpotential, can be labeled as |A^; z/q) c1 , 
where the Ni > and vq E {0, 1} are integers satisfying the constraint 

deg w i 

Y jNl + Y,VQ = N - (i-i^) 

1=1 Q=\ 

The NiS and uqs denote the number of eigenvalues of the matrix X that are equal 
classically to the i th root of W'(z) and to bg, respectively. The U(iV) gauge symmetry 
in a vacuum | Nf, vq) c \ is thus broken down to a product of U(A^) factors. The number 
of non-trivial factors of the unbroken gauge group, i.e. the number of non-zero JVjS, 
is called the rank of the vacuum. In the quantum theory, chiral symmetry breaking 
yields a larger degeneracy of the vacua, that are then labeled as |iVj, fcjj i/q) where 
< ki < Ni — 1. This structure follows from the extremization of the microscopic 
quantum superpotential that will be introduced in the next Section. 
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The expectation values of the operators (1.12) are most conveniently encoded in 
the generating functions 

fc>0 

fe>0 

^hE^^S^- ("«) 

fc>0 

One of our main result is to show that these generating functions satisfy the following 
set of algebraic equations, 

NW' (z)y(z) - y(zf = Aj,(z) , (1.17) 
Nm f ' f ~(z)& f ~ f (z) - y(z)5 f f ' = A f f '(z) , (1.18) 

m f Az)m s ~(z) - y(z)5? = KOz) , (1.19) 

t"'{z)y{z) + NW'{z)M{z) + Nm' f f ~(z)& f ~ f (z) - 2y{z)M{z) = A M (z) , (1.20) 

st ~ ft 

where A y, , A^ and A^ are polynomials. These equations are the famous gener- 
alized Konishi anomaly equations [5], adapted to the case where flavors are included 
in the model [6] and suitably generalized to the extended theory corresponding to 
having arbitrary couplings t. These equations are at the heart of the Dijkgraaf-Vafa 
matrix model formalism [7, 8] , where they follow directly from the planar loop equa- 
tions of the matrix model. They were understood at the perturbative level (i.e. in 
a fixed classical background gauge field) in [5, 6], but a full non-perturbative proof 
requires much more work as explained in great details in [2, 3]. This is where the 
microscopic formalism shows its full power. It is remarkable to reproduce the planar 
matrix model result (which, when flavors are present, also include disk diagrams [8]) 
from finite N gauge theory path integral calculations. As we shall see, these integrals 
can be reduced to non-trivial sums over colored partitions. 

The constraints (1.17), (1.18), (1.19) and (1.20) do not fix completely the expecta- 
tion values. There remains undetermined coefficients in the polynomials that appear 
in the right hand side of these equations. This ambiguity is completely removed by 
the fact that the number of colors iV in the gauge theory is finite and thus only a 
finite number of the operators (1.12) are algebraically independent [9, 10, 11]. Math- 
ematically, this is translated into quantization conditions for the periods of M dz, 

j M{z)Az e2niZ. (1.21) 
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In our microscopic formalism, this is satisfied by construction, since the operators are 
built explicitly from finite N x N matrices and the relations (1.21) will be easy to 
check. 1 

The plan of the paper is as follows. In Section 2, we explain the microscopic 
formalism in the case of the model (1.1). In Sections 3, 4 and 5 we compute the 
scalar, glueball and meson operators respectively. Finally, in Section 6, we extremize 
the microscopic superpotential, derive the anomaly equations and discuss some gen- 
eral features of the solution. We have also included an Appendix containing some 
technicalities used in the main text. 



2 The microscopic formalism 

The starting point of the microscopic formalism [1] is to consider the expectation 
values of the operators (1.12) with fixed boundary conditions at infinity for the adjoint 
scalar field X, 

= diag (a 1 ,...,a N ). (2.1) 

The eigenvalues at infinity 

a = (a 1 ,. . . ,a N ) (2.2) 

are arbitrary fixed complex numbers. We could also try to impose arbitrary boundary 
conditions at infinity for the quark fields Qf and Q* , but we prefer in the present 
paper to first integrate over these fields exactly in the path integrals. The model 
is then reduced to the case with no flavor, but with extra determinant-like factors 
due to the integration over the quarks. The expectation value of an arbitrary chiral 
operator with the boundary conditions (2.1) will be denoted by (a\0\a) and the 
corresponding generating functions in our model are given by 

(a\TrX k \a} 



r^/ x \ a \ IrA a) , . 



fc>0 

fc>0 

fl / )( ,, ) = E (»W. (, 5) 

fc>0 

Clearly, the expectation values {a\0\a) are not physical and the corresponding gen- 
erating functions (2.3), (2.4) and (2.5) do not coincide with the physical generating 



1 In the Dijkgraaf-Vafa matrix model formalism, on the other hand, the conditions follow from the 
extremization of the postulated glueball superpotential with all the required Vcncziano-Yankiclowicz 
terms included [12] and are thus highly non-trivial. 
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functions (1.14), (1.15) and (1.16). 2 In particular, the functions R, S and Gf depend 
on the arbitrary boundary conditions a (as well of course on the other parameters in 
the model), whereas the physical correlators depend on a choice of vacuum but not 
on a. 

The interest in considering the correlators for fixed boundary conditions at in- 
finity is that, at least in an open set in a-space, they can always be computed by 
summing a convergent instanton series [1]. The functions at arbitrary a (outside the 
radius of convergence of the instanton series) are then obtained uniquely by analytic 
continuations. 

There exists a quantum superpotential W m i c (a) for the boundary conditions a [1]. 
One of the fundamental property of this quantum superpotential is that the solutions 
a = a* of the equations 

^f(a = «')=0 (2.6) 

are in one-to-one correspondence with the full set of quantum vacua of the theory 
[1]. The physical correlators in a given vacuum are then obtained by plugging the 
corresponding solution a = a* to (2.6) into the generating functions (2.3), (2.4) and 
(2.5), 

M{z) = R{z;a*) , (2.7) 
Y(z) = S(z;a*), (2.8) 

& f ~ f (z) = G f f (z;a*). (2.9) 

The superpotential W m i c is always unambiguously determined in terms of the expec- 
tation values by a U(1)r symmetry of the model. In our case, the charges for the 
relevant U(1)r are given by 

9 W a X Q f Qf a g m U b t q W (2 1Q) 

U(1) R 11 01 102000002, { ' ) 

where 9 is a superspace coordinate and W an arbitrary superpotential. The Ward 
identity associated with this U(1) R reads 

W mic (a) = J29k d ^- (2.H) 

Using the standard supersymmetric Ward identity 

a\ TrX fc+1 |a) — (k + l) 9 ^ , (2.12) 



2 We have changed slightly our notation with respect to references [1, 2, 3]; R and S in the present 
paper correspond to i? m i c and S mlc in our previous works. 
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Figure 1: The hyperelliptic Riemann surface with the cycles on and chains $ used 
in the main text. Note that the contours a$ are chosen such that they do not encircle 
the points 6q (this is needed in Section 5). The point /iq is taken to infinity. 



we thus get the fundamental formula 

W mic (a) = (a\ TrW(X)\a) (2.13) 

relating the quantum superpotential to the correlators of the chiral operators TrX k . 
We shall also need two additional supersymmetric Ward identities, similar to (2.12), 
that read 

1 (a\TrW a W a X k \a) = N , ( 2 .14) 



16tt 2 x 1 1 1 k + 1 dt k+l 



a\ T QfX k Q f \a) = ^p- (2.15) 

dm, s 



3 The shape function and the scalar operators 

All the chiral correlators for given boundary conditions a can be computed using 
Nekrasov's instanton technology. We plan to provide a general discussion of this 
technology in a purely Af = 1 context in a subsequent paper, but for our present 
purposes all the relevant formulas can be obtained from the existing literature [4, 13, 
14] through rather simple generalizations. 

It is convenient to present first the solution when all the parameters in the problem 
are real, the c^s being widely separated and the 6qs sufficiently smaller than all the 
djS. In this case all the expectation values can be expressed in terms of a so-called 
shape function /, 

{a\t?\a) = 0[f\. (3.1) 

For example, 

R{z-a) = \ [dy^-- (3.2) 
2 Jr z-y 
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The shape function / extremizes the functional 



F=-\J dxdy (x - y) 2 {ln \x-y\- 3/2)f"(x)f"(y) 

i 

+ \JdxJ2(*- b Q ?Hx - b Q ) - 3/2)f"(x) + \ J dxt(x)f(x) (3.3) 



with the constraints 

N 



support[/"] = |J Ii, (3.4) 

£dxf"(x) = 2, (3.5) 
/ dxxf'(x) = 2a; , (3.6) 

N N 

f ( x) = J2\x-a t \, Vx^lJl,, (3.7) 

i=i i=i 

where the N intervals 

h=[w-,wt] (3.8) 

are disjoint. The shape function plays a role that is very similar to the density 
of eigenvalues in planar matrix models. In the present context, the integrals over 
the instanton moduli space are reduced by localization to sums over particular field 
configurations labeled by colored partitions, and the shape function describes the 
dominating colored partition in the limit of vanishing ^-background [13]. A few more 
details are given in the Appendix A. 2. The correlator [a\&\a) for arbitrary complex 
values of the parameters are obtained by analytic continuation, as will be clear in the 
following. 

The formula (3.3) is a simple generalization of cases that were previously studied 
in the literature. The flavorless case of our model was studied for g = in [14] and 
for arbitrary g in [3]. The t = g = theory was studied in [4, 13] in the case where 
the mass matrix rrJ ~{z) is a linear function of z. 

For our purposes, we need to solve the extremization problem in the general case 
where all the couplings are turned on simultaneously. The equation 8J-/8 f(x) = 
reads in this case 



/ dy In \x - y\ f"(y) = £ \n(x - b Q ) + t"(x) , G |J h . (3.9) 
J Q=i i=i 
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Instead of working with /, it is more convenient to study R. Let us note that (3.2) 
implies that / can be obtained from the discontinuity of R across the intervals /«, 

R(x + iO) - R(x - iO) = -inf"(x) . (3. 10) 

The first derivative of (3.9) yields 

1 1 

R(x + iO) + R(x -i0) = ^2 - — — + t"'(x) (3.11) 

Q=l x ~ Q 



U(x) 



+ t"'(x), \/xe\Jl t . (3.12) 



i=i 

This constraint implies that R is a well-defined meromorphic function on the hyper- 
elliptic curve 

N 

V: y 2 = \[{z-wi){z-wr). (3.13) 

i=l 

This curve, with some useful contours, is depicted in Figure 1. In particular, R(z) is 
a two-valued function. From now on, we shall denote by R(z) its value on the first 
sheet, where the asymptotic condition at infinity 

N 

R( z ) = - + 0(l/z 2 ) (3.14) 

z 

is valid, and by R(z) its value on the second sheet. More generally, a hat on a function 
defined on (3.13) will always mean that we consider its value on the second sheet. 
For example, y = —y. Equation (3.12) is equivalent to 

R(z) + R(z) = ^+lf(z). (3.15) 
As any other meromorphic function on the curve R can be written in the form 

R ( z ) = ri ( z ) + r A?l, (3.16) 

where r\ and r 2 are rational functions. Equation (3.15) implies that 

r ^-l{im +t "' {z) )- (3 - 17) 

Moreover, from the integral representation (3.2) and the constraints on the shape 
function /, R cannot have poles on the first sheet. From (3.15), we deduce that the 
only poles of R are on the second sheet at z = bq with residue one. This implies that 
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for some polynomial p(z). The asymptotic behaviour (3.14) implies that degp = 
deg t'" + N + £ and imposes deg t'" + 2 constraints on the coefficients of p. Matching 
the residues at the poles on the second sheet implies £ additional constraints 

p{b Q ) = -\u\b Q )y{b Q ). (3.19) 

There remains N — 1 unknown coefficients in p, as well as the 2N unknown branching 
points wf on the curve (3.13). These 3iV — 1 parameters are fixed by the following 
3iV — 1 independent constraints on the periods of Rdz, 

j> R{z)dz = 2in, (3.20) 

zR{z) dz = 2iirai , (3.21) 

/ R(z) dz = In ( J + NV{no) + 2inZ . (3.22) 

The contours a, and are depicted in Figure 1. The cut-off /i is always understood 
to be taken to infinity at the end of the calculations. Let us note that of all the 
constraints that determine R, none depends on g and thus R itself will not depend 
on g (but will of course depend on a, t and b). Note that the physical generating 
function & will depend non-trivially on g through the solutions a* of (2.6). 

Equations (3.20) and (3.21) directly follow from integrating (3.10) and x times 
(3.10) over the intervals JjS and then using (3.5) and (3.6). Equation (3.22) is more 
interesting. It comes from the integrated form (3.9) of the variational equation for 
the shape function. To see this, let us introduce 



/ R(z')dz' + NlnfiQ. (3.23) 



Since we do not specify the contour used to go from the point at infinity on the first 
sheet /jl to z in (3.23), (3.20) shows that 4>(z) is defined modulo 2i% on the first sheet. 
What happens if we cross one of the branch cuts? Since 0' = R, we can integrate 
(3.12) to find the discontinuity of 0, 

i 

<f)(x + i0) + (j)(x -i0) = Y^ H x - h o) + t"(x) +Ci, Vxeli. (3.24) 

Q=l 

A priori, the integration constants q could depend on the cut Jj. However, by com- 
paring with (3.9), we find that all the qs are actually zero. This means that, by 
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crossing any of the cuts, we go to the same sheet of the function <p. In other words, 
modulo 2in, is well-defined on the curve (3.13). In particular, 

£ 

4>(z) + 4>(z) = H z - h Q) + t"( z ) + 2 ™Z . (3.25) 

Q=l 

This yields 

/ R(z) dz = 0^) - 0(A*o) + 2ivrZ (3.26) 
= -(2N-e)\nfio + t"(fMi) + 2inZ, (3.27) 
which is equivalent to (3.22) thanks to the relation 

t"(z) = NV{z) + \n(U q) , (3.28) 

see (1.4) and (1.8). 

Equations (3.20) and (3.22) imply that 

R(z) dz G 2mZ , (3.29) 

where the integral is computed along any closed contour on the curve ^ (3.13). In 
particular, the function 

F(z) = e^ z) (3.30) 
is well-defined on ^ . Its value on the second sheet is determined by (3.25) to be 

P( z ) _ J(z) _ U.Q=l( Z ~b Q ) t» W _ qU(z) NV (z) (n 01 \ 

The function F has an essential singularity at infinity on the second sheet for 
non-zero V . In the special case V(z) = 0, which corresponds to the conventional 
theory with standard gauge kinetic term, this singularity becomes power-like and F 
is a meromorphic function on . The solution can then be described more explicitly. 
For example, for I < 2N, (3.31) implies that 

F(z) + ^ = H N (z) , (3.32) 

where H N {z) = (1 + qUoS^N) z N + . . . is a degree N polynomial. Equivalently, 

F(z) = l - (h n {z) + y/H N {z)*-4qU{zj) . (3.33) 

11 



Comparing with (3.13), we can relate Hn to the branching points w, 



± 

i > 



N 



H N (z) 2 - 4qU(z) = (1 - qU 5 e , 2N f ]J(z " ™t)(z ~ ™f ) ■ ( 3 - 34 ) 

i=i 

The generating function R then takes the form, 

^ f(z) 2 v(z) yfw^zr - iqu(z) v M ' 2t/(2) ; 1 

To finish this Section, let us comment on the analytic structure of the solution. The 
structure that we have described above is valid for generic values of the parameters, 
but interesting phenomena occur when the boundary eigenvalues of X are chosen 
to coincide with the parameters 6q. By carefully analysing our solution, it is not 
too difficult to show that when 6q approaches a^, the cut Ii = [w~,w^~] closes. At 
a i — i>Q, the curve ^ degenerates to a genus N — 2 curve and the pole at z = bq 
is on the first sheet. More generally, if p distinct a^s are equal to p distinct 6qs, the 
curve degenerates down to genus N — 1 — p and the generating function R then has 
p poles on the first sheet and N — p poles on the second sheet. We have illustrated 
this mechanism on a very simple example in the Appendix A.l. The cases = 6q 
can actually be treated directly and most easily at the level of the sums over colored 
partitions. This is explained in the Appendix A. 2. 



4 The glueball operators 

The inclusion of flavors modifies only slightly the computation of the generating 
function S(z; a) and thus we can follow closely [3]. The fundamental formula relates 
S"(z) to R(z), 

sVH^aM. ( 4 ,) 

fc>l K 

This is the same as equation (3.21) in [3] and the derivation given in that reference 
applies without change when flavors are included. We can also follow closely [3] to 
derive the consequences of (4.1). The only potential difference in the analysis could 
come from the fact that R has poles. However, the residue of these poles are t- 
independent and thus they do not enter in (4.1). From [3] we thus know that (4.1) 
implies that S'(z) must be a meromorphic function on the curve (3.13) of the form 

S'(z) = ^W"(z) + ^ , (4.2) 
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for a certain polynomial s of degree deg W" + iV. The asymptotic condition at infinity 
on the first sheet 

S\z) = O (1/z 2 ) (4.3) 
yields deg W" + 2 conditions on s. Moreover, one must have 



I S'(z) dz = (4.4) 



which yields N — 1 new independent constraints that determine completely s and 
thus S. 

5 The generalized meson operators 

We are now going to show that the generating function for the generalized meson 
operators (2.5) is given in terms of the generating function for the glueball operators 
that we have just computed by the formula 

G f f (z) = ^S(z) nzn f f -^^^Tes w=bQ (m-\w)Y f , (5.1) 

0=1 Q 

where m(z) is the mass matrix polynomial (1.6). Note that all the poles of G f(z) 
are on the second sheet. 

To do the calculation, it is very convenient to use the variations of the functional 
T defined in (3.3). Since T is stationary with respect to the changes of the shape 
function /, we have the simple formula 

bT = — I 54){z)R{z) dz , (5.2) 

where 

i 



W = \H^ Z - H* ~ M " 3 / 2 ) + *(*) • ( 5 - 3 ) 

Q=l 

To derive (5.2), we have used (3.10) and we have defined a to be the sum of the 
contours that circle around the branch cuts of the curve (3.13), 

N 

a = ai . (5.4) 
i=i 

In particular, using (3.21) and (5.2) for 5 = d/dtk, the quantum superpotential (2.13) 
can be rewritten 

fe>0 Ja i=l fe>l ft 



13 



We now use the relation (2.15) combined with (5.5) in the definition (2.5) to obtain 

k>ok'>i z dt k' dm k j 

The partial derivative of T with respect to ^ is then evaluated using (5.2). This 
yields 

nS t \ ST ST 9k' d l dw dip(w) 

G W = LL^f 2^^r~ R{w) ■ (5 - 7) 

k>0 k'>l K Ja kj 

The function ip defined in (5.3) depends on the m{ f only through the bns and t\, 

r k, J 

see (1.9). It is clear that d4>(w)/dm J f does not depend on t k > and thus the partial 

k, J 

derivative with respect to t k > in (5.7) acts only on R(w). We can then rearrange nicely 
the formula using (4.1), 

G f fW = I E 4l / ^^S>) • (5-8) 

The constraint (4.4) ensures that S'(iu) is single-valued along the contour a. We can 
thus integrate by part twice in (5.8) to finally get 

G f M = 4 E "4 / P-^^P-Siw) . (5.9) 

N k z + ^ 2in dm lj 

Now, from (5.3) and using (1.7) and (1.9) we get 

W = 9w^) ; r (510) 

am - am, - 1 

k,f k,f 

The series over k in (5.9) can then be summed up easily and we find 

The resulting contour integral can be computed by deforming the contour a to infinity 
and picking the contributions from the poles at w = z and w = 6q. This yields the 
formula (5.1). 



6 Going on-shell and the anomaly equations 

In this last Section, we are going to solve the equations (2.6). Using (2.7), (2.8) and 
(2.9), we shall then be able to make the link between the off-shell generating functions 
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(2.3), (2.4) and (2.5) that we have computed previously and the physical generating 
functions (1.14), (1.15) and (1.16). In particular, we are going to show that the latter 
satisfy the anomaly equations (1.17)— (1.20). 

The starting point is the fundamental formula that relates the derivative of the 
quantum superpotential to the contour integrals of S'dz, 

dW ■ 1 r 

-^---JfWto + WW-O. (6.1) 

This relation takes exactly the same form as in the theory with no flavor, equation 
(3.51) of [3]. The derivation given in this latter reference, which uses in particular the 
Riemann bilinear relations, applies without modification to the present case. This 
perfect analogy is due to the fact that the poles, that are a priori present in the case 
with flavors, are eliminated when one takes derivatives with respect to a*. 

Let us first examine the consequences of (6.1) for the glueball operators. For 
arbitrary values of a, we have seen in Section 4 that S'(z;a) was well-defined on 
the curve (3.13). Denoting as usual with a hat the value on the second sheet, we 
deduce from (4.2) that 

S'(z; a) + S'(z; a) = NW"(z) . (6.2) 

Integrating, we get 

S(z;a) + S^(z;a) = N(W(z) -W'(pi )) + f S'(z)dz, (6.3) 

where S^(z; a) denotes the analytic continuation of S through the cut Jj. For general 
values of a, S(z; a) is not defined on ^, since the analytic continuation through a 
branch cut depends on the particular branch cut that we choose. However, for the 
particular on-shell values a = a*, the relation (6.1) is satisfied and thus the right- 
hand side of (6.3) no longer depends on i. The physical generating function (2.8) is 
thus well-defined on ^ , with 

y(z)+y(z) =NW'(z). (6.4) 

A trivial calculation using this relation immediately implies that the combination 
A(z) = NW'(z)y(z) — y(z) 2 has no branch cuts, i.e. A(z) = A(z). It cannot have 
poles from the discussion of Section 4. Using the asymptotic condition S(z) = 0(1/ z) 
at infinity, we conclude that it must be a polynomial. This implies the first anomaly 
equation (1.17). 

We can proceed in exactly the same way to derive the other anomaly equations. 
It is straightforward to check that the left hand sides in (1.18), (1.19) and (1.20) have 
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no branch cuts by using (6.4), (3.15) and (5.1) (these last two equations are valid for 
any a, and thus in particular for a = a*). It is also straightforward to check that 
the residues of the possible poles all cancel by using the simple pole structure of the 
various generating functions that we have discussed in the previous Sections. The 
asymptotics at infinity then implies that the right hand sides of (1.18), (1.19) and 
(1.20) must be polynomials. 

Let us close this Section with two remarks. First, we note that the anomaly 
polynomials in (1.17) and (1.20) can be obtained by acting on W mic with first order 
differential operators J n and L n defined exactly as in equations (3.63) and (3.62) of ref. 
[3]. In particular, the Riemann bilinear relations used in [3] to make the derivations 
can be easily generalized to take into account the poles that appear in the generating 
functions R and when flavors are present. On the other hand, to obtain (1.18) 
and (1.19) from variations of the microscopic quantum superpotential, one would 
have to include arbitrary boundary conditions for the quarks in the formalism and 
compute W m i C as a function of both a and these quark boundary conditions. The 
anomaly polynomials in (1.18) and (1.19) would then follow by acting on W m i C with 
suitable first order differential operators containing partial derivatives with respect to 
the quark boundary conditions. In the present paper, we have preferred to integrate 
out the quarks exactly first and thus work with a microscopic superpotential that 
depends on a only. 

Our second remark concerns the set of solutions to the quantum equations of 
motion (2.6). We have shown that any solution must satisfy the anomaly equations 
on top of (3.22) which is valid off-shell. Conversely, the set of solutions to the anomaly 
equations that also satisfy (3.22) is known to be in one-to-one correspondence with 
the full set of quantum vacua of the theory (see for example [10, 11, 12] and references 
therein). One can show that all these solutions also automatically solve (2.6), with 
one rather trivial exception that is discussed below. A simple way to understand this 
point is as follows. First, a straightforward generalization of the analysis in [1] shows 
that vacua of any rank r > 1 of the type \Ni,ki]VQ = 0) are automatically included 
in the set of solutions. Second, one uses the fact that all the other vacua at the same 
rank can be obtained by analytic continuations [11, 12] and thus necessarily solve 
(2.6) as well. 

There is an interesting point concerning the vacua having vq ^ 0. At the classical 
level, one has v q — 1 when one of the a, is equal to bq. At the quantum level, one 
might expect that the solutions to (2.6) associated with these vacua correspond also 
to having = bq. This would be natural from the analysis in the Appendix, that 
shows that if one imposes the boundary condition d = bq, then the quantum function 
R(z; a) has a pole at z = bq on the first sheet. However, what really happens depends 



16 



on the cases one considers and can be more subtle. The subtlety comes from the fact 
that the variables a can undergo non-trivial monodromies, as is well-known from the 
study of the moduli space in the M = 2 supersymmetric theories [15]. Due to these 
monodromies, the actual solution a = a* to (2.6) corresponding to a vacuum with 
vq — \ can actually have all the a^s different from the 6qs. 

The above discussion doesn't apply for the vacua of rank zero. These vacua have a 
completely broken gauge group and correspond to the cases where all the eigenvalues 
djS are equal to the bqs classically. From the discussion in the Appendix, we know 
that the solution is trivial in these cases: the chiral operator expectation values do 
not get any quantum correction. Now, it turns out that these trivial solutions do not 
satisfy (2.6). The reason is that the procedure of integrating out the quarks become 
singular from the point of view of the microscopic quantum superpotential in these 
particular vacua. This can be easily illustrated since these vacua are purely classical. 
Integrating out the quarks from the tree-level superpotential 

W tTee = Tr W(X) + T Qfm f f (X)Q f (6.5) 

amounts to imposing the conditions 

m f f {X)Q f = = T Q f 'm f f -{X) . (6.6) 

The resulting effective superpotential, obtained by plugging (6.6) into (6.5), is simply 

N 

W mic = Tr W(X) = W( ai ) , (6.7) 
i=i 

whose variations only yield W'(a,i) = 0. For these solutions, the matrix m is 
invertible and Qj = = 0. The superpotential (6.7) is thus missing the solutions for 
which m j has zero eigenvalues and T Q^Qf ^ 0. These solutions correspond precisely 
to the cases Oj = 6q. When the rank of the solutions is r > 1, and contrary to the 
case r = 0, there are non-trivial quantum corrections and as we have explained above 
the solutions are actually obtained from (2.6). 
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Appendix: the special cases with a$ = 6q 
A.l A simple example: N = Nf = 1 

Let us consider the solution for the generating function R(z; a) described in Section 
3 in the case N — Nf — 1, U = z — b and V = 0. If we note a± — a and bi = b, we get 

m = - f 1 * + 1 = ( i - Z 7 a+ J ) (A.i) 

from (3.35) and (3.21). The function R(z) is two- valued, with asymptotics at infinity 
R(z) ~ l/z on the first sheet and a pole at z = b on the second sheet. The branching 
points are given by the equation 

(z - a + qf - Aq(z - b) = (z - w~)(z - w + ) , (A.2) 

which yields 

w ± = a + q ± 2^q(a - b) . (A.3) 
When b — > a, we see that the two branching points collide and R reduces to 

*M = ^- (A.4) 

As explained in the main text, this is a very general phenomenon: when 6q — > a^, 
the cut closes and the pole at z = b Q is on the first sheet. We are going to 

find this property again in the next subsection from a direct analysis of the sum over 
colored partitions. 

A.2 Direct analysis using the sums over colored partitions 

Chiral correlators with fixed boundary conditions a for the field X are given by a 

— * 

sum over colored partitions k of the form 

(a\e\a) =lim Sffl^ , (A.5) 

where the limit e — > corresponds to a vanishing ^-background. A colored partition 
is a collection k = (ki, . . . , k^r) of ordinary partitions k« = {k^ a }, fc^i > fcj^ > • • • > 
fcj > 0. Much more details about these sums can be found for example in [2]. The 
measure over the set of colored partitions decomposes as 

^ = /l/ijW,eK f (a,e), (A.6) 
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where the dressing factor gives the contribution from the integration over the quark 
fields, 

TV &i,a 

^=nnn f/ (^+ e (' 9 -«))- ^ 

i=l Q =l /3=1 

This formula generalizes the dressing factor obtained in [4] in the case of a linear mass 
function m ^(z). The polynomial U is defined in (1.7). When e — > 0, the sums (A. 5) 
are dominated by a single large colored partition described by the shape function / 
[13]. This has been used extensively in Section 3. 

The cases where some of the OjS are equal to the b Q s are special. What happens 
is very clear from the form of the dressing factor (A. 7): if a, = 6q, only the trivial 
partition kj = yields a non-zero contribution. In general, if p distinct a^s are equal 
to p distinct 6qs, the sum over colored partitions reduces to a sum over N—p ordinary 
partitions which can be computed as in Section 3. The dominant colored partition is 
described by a smooth shape function / that extremizes a functional given by 

F=-\J dxdy (x-y) 2 (\n\x-y\- 3/2) f"(x)f"(y) + ^ J dxt(x)f"(x) 

+ fdx(x-b Q ) 2 {ln(x-b Q )-3/2)f"(x), (A.8) 

Q=l J 

where vq — 1 if &q = aj and vq = otherwise. The constraints on / are similar to 
(3.4)-(3.7), except that now the support of / is made of N — p distinct intervals cor- 
responding to the GtjS that are distinct from the 6qs. One can solve this extremization 
problem as in Section 3. The resulting generating function 

R(z) = y^- + - [dy^l (A.9) 
1 ; ^z-b Q 2j R y z-y y } 

is defined on a hyperelliptic curve of genus — 1 — p, with p poles on the first sheet 
and N — p poles on the second sheet having residue one. As already emphasized, 
this is exactly the same solution as the one obtained starting from R(z; a) for generic 
values of a and b and then going to the special points a« = 6q. 
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